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ABSTRACT 

We analyze the evolution of binary stars to calculate synthetic rates and delay times of the most 
promising Type la Supernovae progenitors. We present and discuss evolutionary scenarios in which 
a white dwarf reaches the Chandrasekhar-mass and potentially explodes in a Type la supernova. We 
consider: Double Degenerate (DDS; merger of two white dwarfs), Single Degenerate (SDS; white dwarf 
accreting from H-rich companion) and AM Canum Venaticorum (AM CVn; white dwarf accreting from 
He-rich companion) scenarios. The results are presented for two different star formation histories; 
burst (elliptical- like galaxies) and continuous (spiral-like galaxies). It is found that delay times for 
the DDS in our standard model (with common envelope efficiency acE = 1) follow a power-law 
distribution. For the SDS we note a wide range of delay times, while AM CVn progenitors produce a 
short burst of SNe la at early times. The DDS median delay time falls between ~ 0.5—1 Gyr; the SDS 
between ~ 2 — 3 Gyr; and the AM CVn between '-^ 0.8 — 0.6 Gyr depending on the assumed qce- For 
a Milky Way-like galaxy we estimate the rates of SNe la arising from different progenitors as: ^ 
yr^i for the SDS and AM CVn, and ~ 10"^ yr"^ for the DDS. We point out that only the rates for 
two merging carbon-oxygen white dwarfs, the only systems found in the DDS, are consistent with the 
observed rates for typical Milky Way-like spirals. We also note that DDS progenitors are the dominant 
population in elliptical galaxies. The fact that the delay time distribution for the DDS follows a power- 
law implies more Type la supernovae (per unit mass) in young rather than in aged populations. 
Our results do not exclude other scenarios, but strongly indicate that the DDS is the dominant 
channel generating SNe la in spiral galaxies, at least in the framework of our adopted evolutionary 
models. Since it is believed that white dwarf mergers cannot produce a thermonuclear explosion given 
the current understanding of accreting white dwarfs, either the evolutionary calculations along with 
accretion physics are incorrect, or the explosion calculations are inaccurate and need to be revisited. 
Subject headings: binaries: close — supernovae: general 



1. INTRODUCTION 

Type la Supernovae (SNe la) play an important role 
in astrophysics as cosmological distance indicators. Ad- 
ditionally, they provide iron peak elements, having di- 
rect consequences f or the ch emical evolution of galaxies 
jRiess et al. 1995: Matteucci & Greggio 1986: de Donder 
fc Van beveren 2003). Currently > 2000 SNe la have been 
observationally confi rmed^, some as distant as z ~ 1.55 
(Strolger et al. 20041. Empirically-derived relationships 



between ligh t curve 
(i.e.. Amis 



proper ties and intrinsic lumi nosity ^"^S! |Thielcmann et al 
HnnQh f^^i.^^ „ /lT3„„i I which wouldTje produc 



that the physical properties of their progenitors are un- 
changing with redshift. 

Despite the continued use of SNe la as standard can- 
dle distance indicators, their origin remains uncertain. 
It is generally accepted that SNe la arise from the total 
disruption of a Chandrasekhar mass 1.4 M©) carbon- 
oxygen white dwarf (WD) as a result of thermonuclear 
explosion. This hypothesis is supported by the fact that 
the amount of energy ob served in the explosions 10^^ 



mutter et al. 



19931 and stretch-factor s, (Perl- 
lave made it possible to 'standard- 
ize' absolute magnitudes of SNe la light curves over a 
wide variety of host galaxy environments.^ Their use as 
'standard candles' on cosmological scales has led to the 
realization that the expansion rate of the universe is ac- 
celerating, an d has enab led accurat e estimations of f^A 
and Q.M (e.g. jSchmidt et al..l998,.Riess et al.[1998||Perl- 
mutter et al.||T§99). However, using SNe la m order to 



is equal to the amount 
which would be produced in the conve rsion of carbon and 
oxygen into iron (see e.g., Livio|2000 for a review). It is 
natural to presume that the exploding white dwarf must 
accrete matter from a close stellar companion until reach- 
ing the critical Chandrasekhar mass, though the nature 
of the companion, the rate and efficiency at which mass 
is accumulated onto the white dwarf, and which array of 
conditions are necessary in order for the whi te dwarf to 
ignite ex plosively, are not well understood (Nomoto et 



set the distance scale for the determination of cosmolog- 
ical quantities requires the (still unfounded) assumption 

^ http:/ /www. cfa.harvard.edu/iau/lists/Supernovae. html 
^ It has been shown that SNe la originating among young stellar 
populations are overall more luminous than those associated with 
older stellar populations (i.e., |Hamuy et al.|1995| l. 



al. 19971 



order to constrain the nature of SN la progenitors, 
the SN la rate has been studied as a function of parent 
galaxy stellar mass, star formation rate, colo ur, morphol- 



ogy, and radio power by several groups (e.g., Mannucci etl 



al.i|2005i iSullivan et al. 


120061 ICalura & Matteucci 2006f 


Delia Valle & Panagia 


2003 


1. Mannucci et al. 


( 


2005k 
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found that the SN la rate is higher in bluer, later type 
galaxies, supporting the hypothesis that there is a non- 
negligibl e number of SNe la or iginating from young pro- 
([2006| found that the SN la rate 



Sullivan et al. 



genitors. 

per unit mass increases as a function of star formation ac- 
tivity a trend which coincides with the results of.Calura, 



fc Matteucci ( 2006 1 , who used chemical evolution mod- 
els to derive t herate of SNe la as a functio n of galaxy 
Hubble type. IDella Valle & Panagial (pOOSl discovered 



an enhanced SN ia rate in radio-loud galaxies compared 
to their radio quiet counterparts, possibly as a result of 
past interactions/mergers with dwarf galaxy companions 
leading to an increased n umber of newly-formed or cap- 



tured young stars (see also Delia Valle et al. ( 2005[ )). The 



evolution of the supernova rate (both core-collapse and 
Ty pe la) as a function of cosmic time was investigated 
by 
ica 



Madau et al. 



( 1998 1, who convolved a set of theoret- 

characteristic SN delay times with the cosmic star 
formation history, resulting in an estimate of SN rates 
out to intermediate redshifts. Later studies showed that 
a single-component delay time could not be reconciled 
with the observed mass of iron in galaxy clusters a nd the 



corresponding rat io of core-collapse to la SN rates ( Maoz 
fc Gal-Yam|[2004l ) 



SNe la are observed in both y oung and old galaxies 
(Branch & van den Bergh 19931, thus it is natural to 
presume that the progenitors may originate from both 
young and old stellar populations. Recently, it has been 
found that SNe la appear to span a wide range of delay 
times which is bimodal in nature, consisting of a 'prompt' 
population with short delay times, and a 'tardy' popu- 
lation whose average delay time distribution (DTD) is 
much wider and is best desc ribed by a decaying exponen- 
tial funct ion ([S cannapicco fc Bildsten|2005 Mannucci et 



al. 2006 



parent bimodal DTD shape is limited to low-redshift SNe 
alone or whether i t also applies to SNe at higher redshift 
is still unclear ( ^Dahlen et al.| 2008|. 

Two formation scenarios nave emerged as the most 
likely channels for SNe la progenit ors: The Single De- 
generate Scenario (SDS, Whelan fc Iben||1973; Nomoto 
1982[) and the Doub le Degenerate Scenario (DDS7 |lben fc 
Tutukov|1984| fVV^bink,, 19841. The SDS is encountered 



Dilday et al. 112008). However, whether the ap- 



when a WD accretes H-rich matter during stable Roche 
Lobe Overflow (RLOF) from a stellar companion; either 
a main sequence or an evolved (giant) star. The WD 
increases in mass up to the Chandrasckhar mass limit, 
enabling carbon to ignite explosively in the WD center 
causing a SN la. The DDS is the result of a merger of 
two white dwarfs. If the combined mass of the merger 
exceeds the Chandrasckhar mass, the result may be a SN 
la. Additionally it has been suggested that a third chan- 
nel, the AM CVn channel'^, m ay account for 1 % of SNe la 
dSolheim fc Ytingelson|2005t. AM CVn binaries (see e.g., 
|Nelemans |2005[ I Warner|1995p are ultra-compact systems 
involving a WD accretor and a helium-rich donor ex- 
changing matter via RLOF. The donors are expected 
to be small stars given the small orbital size (close or- 
bits; < 1 hr). Other possible SNe la formation scenarios 
have been proposed, though they hkely do not account 
for the majority of SNe la (e.g., s ub-Chandrasekhar mass 
SNe la ( Woosley fc Weaver|1994 1 , common envelope WD 



mergers (|Livio fc Riess||2003 Applegate 1991 Sparks fc 
|Stecher|[lF7ll) ]: 

In this study we follow the evolution of stellar popula- 
tions in two different environments, that which is typical 
for an elliptical galaxy (instantaneous burst of star for- 
mation at t — 0) and that typical for a spiral galaxy 
(continuous star formation). We also employ two differ- 
ent parameterizations for common envelope (CE) evolu- 
tion. To substantiate our conclusions, we use exactly the 
same total stellar mass and metallicity for each popula- 
tion; the only differences in our model galaxies are the 
star formation histories, and the assumed CE removal 
efficiency acE- We show which SN la progenitors are 
the most likely (from an evolutionary perspective) for 
each host galaxy type, and derive delay times and rates 
for the most promising SNe la progenitor scenarios. We 
compare our results to those of previous studies and we 
discuss our results (e.g., in terms of explosion physics) in 
the last two sections, respectively. 

2. MODEL DESCRIPTION 

Our stellar evolution calculations are performed using 
the Star Track population synthesis code. A detailed de- 
script i on of the in put physics is presented in |BelczynsE] 
et al. (2002 2008). Single star evolution is followed from 



the !ZAMS until remnant formation employing modified 
analy t ic for mulae and evolutionary tracks from |Hurley| 



et al. (20001). Evolution of binary stars is more complex, 
and several processes important for field binary evolu- 
tion are accounted for, such as tidal interactions, mass 
transfer phases, common envelope evolution, supernova 
kicks, magnetic brak ing and gravitational radiation (see 
Belczynski et al.||2008, for formulae). We incorporate re- 



cent prescriptions tor mass growth of WDs, employing 
accretor mass-dependent accumulation efficiencies which 



thick WD winds (Nomoto et al.I'2007; Kato fc Hachisu 
2004.,Hachisu et al.,1999||Kato fc Hachisu 1999; PrialnIC 


fc Kovetz 


1995 Hashimoto et al. 


1986 see below). The 



the evolution. 

A merger between two white dwarfs may lead to a DDS 
SN la. A detached WD-WD binary will eventually reach 
contact due to angular momentum loss from the emission 
of gravitational radiation, and if the binary configuration 
(e.g., mass ratio) leads to a merger, the less ma ssive WD 
is accreted o nto the more massive WD (see e.g., Tutukov 
[fc Yungelson||1979 ) . If the combined mass of the merger 



exceeds the Chandrasckhar mass, it is recorded as a po- 
tential DDS SN la progenitor. We assume a priori that 
every WD-WD merger with M > 1.4 Mq consisting of 
CO-CO, CO-He, or He-He WDs leads to an instanta- 
neous SN la.^ Mergers involving other WD types occur 
in our simulations but we do not count them as SN la 
progenitors. 

A SDS or AM CVn SN la may result from the accu- 
mulation of matter on a white dwarf's surface via sta- 
ble RLOF from a stellar companion. For accretion of 
hydrogen-rich material, strong nova explosions inhibit 
the accumulation of hydrogen on the WD surface for 



* In 

tween 



SWB-Iike in Belczynski et al. ( 2005 I 



[Belczynski et H] |2005t , it was assumed that a merger be- 
any two WDs, mcludmg ONe WDs, witli a combined mass 
exceeding 1.4 Mq would lead to a SN la, though CO-CO mergers 
made up 88 % of the DDS SNe la in that study. 
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very low ma ss transfer rates < 10 yr (Prialnik &; 
Kovetz|[l995| . For mass transfer rates above this thresh- 
old, w e mterpolate over the results of |Prialnik fc Kovetz 



(19951 to obtain the mass accretion efficiencies, and we 



accou nt for optically thick winds at high mass transfer 
rates (Hachisu et al. 1999, see also Belczynski, Bulik & 
Ruiter (2005) Section 2.3 for a more detailed descrip- 
tion). The only difference between the accumulation 
efficien cies in this work and those of [Belczynski et atL\ 
(|2005) is that here we have additionally included an up- 
dated piescri£tion_for accretion of hydrogen-rich matter 



from pomoto et al. (2007), in which fully efficient ac- 
cumuiation is only achieved for a very narrow range of 
mass transfer rates, a nd is also dependent u pon the white 
dwarf accretor mass (Nomoto et al. 2007 see equations 
5 & 6). 

Accretion of helium-rich matter is treated in the same 
fashion as in Belczynski et al .l (|2005l) : accretion prescrip- 
tions are adopted from Kato & Hachisu (1999). However 



this study contains one major difference which affects our 
results: though in this work we do allow for the formation 
of sub-Chandrasekhar mass SNe la, we do not include 
these binaries as pot ential SN la progenitors here. In |Bel-| 
czynski et al. (2005), a large fraction (61 %) of the Type 
la SMe which contributed to the presented delay times in 
the standard model were in fact sub-Chandrasekhar SNe 
la, in which the accumulation of ~ 0.1 Mq of He-rich 
material on the WD surface could lead to an edge- lit 



detonation and subsequent SN la (Woosley & Weaver 
19941 IKato & Hachisul 19991) 



The criteria used here for defining SNe la progenitors 
arising from different formation channe ls are different 

In this work, 



(2005) 



from the work of |Belczynski et al. 
we only consider accreting WUs which have obtained a 
mass of 1.4 M0 as potential SNe la. We make the dis- 
tinction between SN la progenitors with CO white dwarfs 
accreting from non-degenerate hydrogen-rich compan- 
ions (SDS), and white dwarfs accreting from helium-rich 
companions (AM CVn). We note that He WDs never 
reach the Chandrasekhar mass in our simulations with 
the adopted accumulation physics, and we have assumed 
that oxygen-neon-magnesium WDs collapse to form a 
neutron star (accretion induced collapse) upon reach- 
ing the Cha ndrasekhar mass, r ather than producing a 



SN la (e.g., [Miyaji et al.|1980[ ). For the SDS and AM 

CVn cases, we record the binary as a SN la once the ac- 
creting WD has reached the Chandrasekhar mass. The 
SDS may occur when a WD accretes matter via RLOF 
from any hydrogen-rich companion (e.g., main sequence 
or evolved star) . If the WD accumulates enough hydro- 
gen on its surface such that steady burning can occur, the 
WD can increase in mass up to 1.4 Mq, carbon ignites 
explosively in the WD centre and the result is a SN la. 
In the AM CVn scenario, we assume the result is a SN 
la if the CO WD reaches 1.4 M© via stable RLOF from 
a helium-rich companion. We allow for different types 
of helium-rich donors in the AM CVn scenario: helium 
stars (non-degenerate stars burning helium in the core 
or in a shell which have been stripped of their outer hy- 
drogen envelope), helium white dwarfs, and hybrid white 
dwarfs (CO-rich core, helium-rich mantle). 

We adopt two contrasting star formation rates (SFR): 
an instantaneous burst at t = (elliptical galaxy) and a 
constant SFR for 10 Gyr (spiral galaxy). Both popula- 



tions are then evolved up to 15 Gyr. The mass formed 
in stars in both cases is the same: 6 x 10^° M©, which 
corresponds to the stellar mass in the Milky Way (MW; 
Klypin et al. 2002 1. In each population we adopt a binary 



fraction of 50% (2/3 stars in binaries), though this frac- 
tion may be overestimating the binary population among 
low-mass stars ( Lada|2006 1, and un derestimating the bi- 
nary fraction among massive stars (Kobulnicky & Fryer 
2007). All stars are evolved with solar-like metallicity 
{Z = 0.02). While the initial distributions represen- 
tative of the physical characteristics of ZAMS binaries, 
and the correct way in which to treat common envelope 
evolution and magnetic braking are all somewhat uncer- 
tain, our choices for various distribution functions are 
constrained by available observations (see below). The 
m agnetic braking pr e script ion which we adopt is that 
of [Ivanova fc Taam (2003), which is based on a two- 
component coronal model (see Section 3.2 of Belczynski 
jet al. 2008). 

ZAMS masses (Mzams) span the entire mass range: 
0.08 — 150 M0. Single stars and binary primaries are 
drawn fro m a 3-component br oken power law initial mass 
function ( jKroupa et al.|1993 l, and secondary masses are 
obtained from a hat mass ratio distribution (Mazeh et 



all (1992 1, q ^secondary /primary), which is th e canoni- 
cal ch oice among population synthesis studies (Lii et al. 
2006|. However, given the observational selection ettects. 



the true mass ratio distribution among ZAMS binaries 
remains unk nown, though it is likely dependent upon 
stellar mass ( |Trimble|[l990| |Duquennoy fc Mayor''1991). 
It has been suggested that the mass ratio distribution 
among local spectroscopic field binari es as well as young 



early-type stars is peaked near unity ( Fisher et al 



2005 



Kobulnicky fc Fryer||2007" I . Initial orbital separations m 



our calculations span a wide range up to 10 R© and are 
drawn from a distribution oc l/a ( Abt|| 198*3 ), which has 
been found to be rep resentative of the local pop ulation 
of Hipparcos binaries ( Lepine & Bongiorno |2007l) . Initial 
eccentricities are drawn fro m a distribution ^(e) = 2e 
(Duq uennoy fc Mayor||1991[ ). 

Close binaries (and potential SN la progenitors) are 
believed to encounter a phase of common envelope evolu- 
tion, and this remains one of the most poorly-understood 
phases in stellar astrophysics. For this reason, we present 
the delay times and rates of SNe la for two different CE 
parameterizations. Some comparison between different 
prescriptions for CE evolution have been t ested against 
observation s for local double white dwarfs (Nelemans & 
Tout 2005 ) . It is unclear at this point how the com- 
mon envelope phase should best be treated in population 
synthesis studies, and currently detailed models arc not 
sophistica ted enough to explore the parameter space in 
detail (Ricker fc Taam 2008 j ). For th is work, we cho ose 
the 'a prescription for CE evolution ( Webbink|1984 1, in 
which the orbital energy of the binary is used to unbind 
the common envelope from the system. We choose two 
different values for the parameterization of common en- 
velope removal efficiency acE- For our standard model 
(Model 1) we choose acE x A = 1, where A is a function of 
the donor envelope structure, a nd is of order unity (see 
also van der Sluys et al. |2006 1. As an alternative, we 
additionally include a model with decreased common en- 
velope removal efficiency in which acE x A = 0.5 (Model 
2). 
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3. DELAY TIMES 

Model 1: acE x A = 1. We use the elliptical model, 
with all stars born at t = to demonstrate the delay time 
distribution for the various progenitors. Figure 1 (top 
panel) shows the characteristic delay times for Model 1, 
and the average and median delay times of all three pro- 
genitors are also indicated. The sharp cut-off near 15 Gyr 
is artificial, as evolution was only allowed to proceed for 
15 Gyr. 

The DDS events follow a power-law like DTD, with a 
median of iMed = 0.93 Gyr and approximated functional 
form of f {t) (X 100t~^. Only a small fraction ~ 5% are 
'prompt' (Mannucci ct al. 2006[) events with delay times 
t < 100 Myr. The DDS events are expected to be found 
long after the star formation has ceased (t ^ 10—15 Gyr) 
therefore we expect these progenitors to produce SNe la 
both in young (spiral and starburst) and old (elliptical 
and bulges) host populations. Even though we have al- 
lowed DDS progenitors to arise from mergers between 
any combination of CO and/or He WDs, all of our DDS 
systems originate from mergers of CO-CO WD binaries, 
since mergers of He- He or CO-He WD binaries never ex- 
ceed 1.4 M0. The evolution leading to the formation of a 
CO-CO WD binary usually starts with two intermediate- 
mass stars (MzAMS ~ 3 — 9 Mq) that evolve through a 
series of close interactions, the first one typically being 
stable RLOF and the last one being a CE phase. Once 
a CO-CO WD is formed, the dominant mechanism for 
angular momentum loss is g ravitational r adiation, the 



merger timescale Tgr oc a (Peters 19641. Given the 



initial distribution of orbital separations cx a~^, and the 
evolutionary orbital change that is to first order the same 
for all DDS progenitors (reduction of orbital size by a fac- 
tor of ^ 10 — 100 during the CE phase), the delay time t 
should likely then follow from the product of the initial 
distribution and the change in separation which occurs 
due to gravitational radiation: 

a-^da/dt cx (ii/^)-!*-^/* oc t"^. Thus, it is not sur- 
prising that the DTD follows a power law as presented 
in Figure 1. 

It is found that relatively speaking, the shortest de- 
lay times originate from the AM CVn channel. The 
distribution is somewhat narrow, with the majority of 
events occurring with t < 2 Gyr, with a median of 
tucd = 0.59 Gyr. Most AM CVn progenitors originate 
from intermediate-mass stars (primaries Mzams ~ 5 — 7 
M0, secondaries Mzams ~ 2 — 4 Mq), however, in con- 
trast to DDS progenitors, they undergo two CE events. 
After two CE phases the orbits of pre- AM CVn binaries 
are very close, which allows for rapid orbital decay due to 
emission of GR, and the final RLOF starts with no sig- 
nificant delay. At the onset of the final RLOF (start of 
AM CVn phase), the accretor is a massive CO WD and 
the companion is either a helium star (65 %) or a he- 
lium/hybrid WD (35 %). Once the RLOF commences, 
it proceeds on a short-timescale with high mass trans- 
fer rates (- 10"^ - 10"^ M© yr^^), enabling the CO 
WD to rapidly build up to the Chandrasekhar mass and 
subsequently produce a SN la explosion. Note that the 
potential SN la progenitor systems discussed here in the 
framework of the AM CVn scenario are by no means a 
representation of the observed sample of ~ 22 AM CVns. 
The rare systems discussed here are more massive, ultra- 



compact systems (median Porb ^ 10 min) and thus have 
higher mass transfer rates than the population of AM 
CVns which is presently observed (median Porb ~ 35 
min). The typical, low-mass AM CVn binaries are long- 
lived systems which do not disappear from the observa- 
tional population due to SN la disruptions. There is a 
clear observational bias against those AM CVn binaries 
which can potentially produce a SN la. 

The SDS channel displays a rather fiat (though some- 
what decreasing) DTD over a wide range t ^ 1 — 15 
Gyr, with a median of tucd — 3.23 Gyr. The SDS sys- 
tems are found to be binaries with a CO WD accreting 
from an evolved star 95%; mostly red giants) or a 
main sequence star ('^ 5%). Since the evolution leading 
to the formation of a CO WD is rather fast, the delay 
time for the SDS is set by the evolutionary timescale 
of the evolved donor. That in turn is a strong func- 
tion of initial ZAMS donor mass (0.7 <Mzams < 2.7 
M0; lifetimes 0.5 — 30 Gyr). The longest delay times fol- 
low from progenitors whose donors, on average, had the 
smallest initial masses. Mass transfer rates can be ini- 
tially high; up to ~ 10^^ Mq yr~^, in the case of giants 
with Mzams > 1 M©, though hydrogen accumulation 
is only fully efficient for a narrow range of mass trans- 
fer rates (see Nomoto et al.|[2007 Figure 4) . Despite the 
high mass transfer rate frorri the donor, only a fraction of 
the hydrogen accumulates on the white dwarf (accumu- 
lation rates of ~ 10~^ Mq yr~^). For SDS progenitors 
with main sequence donors, the two stars are brought 
into contact due to loss of orbital angular momentum 
from magnetic braking (convective secondary) and to a 
lesser extent, gravitational wave emission. Mass transfer 
and accretion rates are much lower (~ 10"^"'^ Mq yr~^), 
and the WD takes > 1 Gyr to accrete up to the Chan- 
drasekhar mass. 

Model 2: acE x A = 0.5. The delay time distribution 
for Model 2 is shown in Figure 1 (bottom panel). The 
change in the treatment of the CE evolution leads to 
smaller orbital separations after the common envelope 
phase, which affects the subsequent binary evolution and 
thus leaves an imprint on the resulting DTD. 

As with Model 1, Model 2 DDS progenitors involve 
only CO-CO WD binaries. However with the decreased 
CE efficiency, many potential DDS progenitors merge be- 
fore a detached white dwarf binary has been produced, 
thus there is a lack of DDS progenitors and there is an 
increase in the number of merging WD + AGB core sys- 
tems. The same power-law curve from the top panel of 
Figure 1 is shown for comparison. The power-law-like 
shape of the DDS is present in Model 2 for delay times 
< 6 Gyr with an additional pile- up of progenitors at short 
delay times (t ^ 1 Gyr). The pile-up is due to the fact 
that in Model 2, DDS progenitor systems which survive 
the CE phase (e.g., binaries which do not merge in the 
common envelope) are on closer orbits upon emerging 
from the CE, thus they merge with relatively shorter de- 
lay times {tgr oc a**). This 'shift' to earlier delay times 
in the DTD from Model 1 to Model 2 serves to build 
a relatively stronger peak at short delay times, while at 
the same time decreases the number of DDS progenitors 
with delay times > 1 Gyr. This is the reason for the 
shorter median delay time in Model 2: ^Mcd = 0.52, vs. 
that of Model 1 (tMcd = 0.93). 

The SDS channel of Model 2 is more efficient than 
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that of Model 1 since the WD and the stellar compan- 
ion are brought on a closer orbit during the CE phase. 
For this reason, wider systems which would have evolved 
to become double white dwarfs (e.g., AM CVn double 
degenerates) in Model 1 evolve into binaries with WDs 
accreting from non-degenerate companions in Model 2. 
In Model 2, the contribution of main sequence donors in 
the SDS channel is slightly higher since after the common 
envelope and formation of a CO WD, in many cases the 
stars are close enough for the secondary to fill its Roche 
lobe before evolving off of main sequence. The number of 
SDS progenitors involving main sequence donors is dou- 
bled relative to Model 1 and now constitutes 10 % of 
the total SDS population (90 % are giant or sub-giant 
companions). SDS progenitors originate from systems 
with donor masses ^ 0.7 < Mzams < 2.7 Mq, many of 
which are Mzams > 1-25 Mq thus have main sequence 
lifetimes shorter than ~ 5 Gyr. The increased number of 
SDS progenitors at earlier delay times leads to a decrease 
in the SDS median delay time from ^ 3.2 Gyr (Model 1) 
to ~ 2.1 (Model 2). 

The DTD of the AM CVn progenitor population of 
Model 2 is notably different from that of Model 1, in 
that it is bimodal. The 'fast' channel at delay times < 3 
Gyr is still present, and originates from progenitors with 
donor masses Mzams }^ 2.5 Mq which undergo two CE 
events. However in Model 2, the decreased orbital sepa- 
ration of post-CE binaries allows for the formation of a 
new SN la progenitor channel, in which a CO WD ini- 
tially accretes hydrogen from a low-mass main sequence 
star for several Gyr, as the main sequence star continues 
to fuse hydrog en into helium in its core (see also Podsi- 
adlowski|2008p . The main sequence star, which has been 
losing its outer layers in RLOF, has built up a signifi- 
cant helium core by the time its mass is depleted to the 
hydrogen burning mass limit (0.08 Mq). Upon reaching 
the hydrogen-burning mass limit, the donor is said to be 
degenerate, and since it is helium-rich (exposed core) we 
treat it as a helium white dwarf. There is a brief period 
(several Myr) where RLOF ceases due to the decreased 
size of the newly degenerate helium WD before the stars 
are brought in contact again due to emission of gravita- 
tional radiation. Stable RLOF begins between the CO 
WD and helium WD, until the CO WD accretes up to the 
Chandrasekhar mass, exploding in a SN la. This addi- 
tional (and scarcely-populated) 'slow' AM CVn channel 
with delay times > 6 Gyr originates from progenitors 
with donor masses 0.7 < Mzams < 1-0 Mq, which only 
undergo one CE event, and evolve through an active cat- 
aclysmic variable (CV) phase for > a few Gyr. Under the 
Model 1 CE evolution, such a system would be found as 
a detached COWD -I- MS binary, and would not make a 
CV, nor a SN la, in a Hubble time. 

4. RATES 

Assuming a binary fraction of 50 %, as we have done in 
this study, 0.17 % and 0.09 % of stellar systems (where 
a 'system' represents either a single star or a binary) 
evolve into SNe la in Model 1 and Model 2, respectively. 
For an assumed binary fraction of 100 %, 0.34 % and 
0.18 % of binaries evolve into SNe la in Model 1 and 
Model 2, respectively. Note that this applies to both 
elliptical and spiral populations as they differ only in 
SFR but not in any evolutionary parameters other than 



acE X A. The above fractions translate into integrated 
efficiencies of 1 SN la per 2500 Mq and 4700 Mq of 
formed stars for Models 1 and 2, respectively, when a 
binary fraction of 50 % is assumed. For an assumed 
binary fraction of 100 % the integrated efficiencies are 1 
SN la per 1500 Mq and 1 SN la per 2800 Mq for Models 
1 and 2, respectively. Though the efficiency scales with 
the adopted binary fraction, the rates as a function of 
time, obviously, are different depending on the adopted 
SFR. 

A general summary of Model 1 and Model 2 rates is 
presented in Tables 1 & 2 where we show the rate of SNe 
la in S NuM, where 1 SNuM = 1 SN (100 yr)-i (10^° 
MfTi)~^ (|Mannucci et al.||2005|). We show the number 



of SNe la in SNuM for the three Chandrasekhar-mass 
models investigated in this work at 4 different epochs: 
0.5, 3, 5 and 10 Gyr. SNuM rates are shown for our 
four galaxy models: elliptical. Model 1; spiral. Model 1; 
eUiptical, Model 2; spiral. Model 2. 

Model 1. In Figure 2 (top panel), SN la rates (number 
of SNe per unit time) are shown for our Model 1 elliptical 
galaxy (acE x A = 1, instantaneous starburst at t = 0). 
The three progenitor types: DDS, SDS and AM CVn 
are shown separately. The DDS rate declines with time, 
but DDS progenitors are found both at early and late 
times. The DDS events dominate over the other progen- 
itor types. This dominance is very strong (^1—2 orders 
of magnitude) and holds for the entire evolution of an 
elliptical galaxy (0 — 15 Gyr). The rate of SNe la from 
AM CVn systems is high for times < 2 Gyr and then 
rapidly declines. Since for a typical AM CVn progenitor 
(i) the CO WD forms early (tcvoi < 100 Myr) and with a 
high mass (Mcowd 1.1 Mq; from primaries of initial 
mass Mzams ~ 5 — 7 M©), and (ii) the orbital separa- 
tion after two common envelopes is small^, the two stars 
are brought into contact either by GR (He WD) or a 
combination of GR and evolutionary expansion (helium 
star). Once the donor fills its Roche Lobe, accumula- 
tion is fully efficient for binaries with helium star donors 
(mass transfer rates ~ 10~^ Mq yr~^), though in the 
case of helium WD donors, the mass transfer rate is ini- 
tially higher (10"^ - 10"^ Mq yr'^) due to the smaller 
stellar separation of these systems upon reaching contact 
0.3 vs. 1.0 Rq). For these relatively high mass trans- 
fer rates, only a fraction (~ 50 %) of the mass is accreted 
by the CO WD. In either case, the mass accretion rate 
averaged over time from contact to SN la is found to be 
on the order of 10"^ - 10"^ Mq yr"! for AM CVn SN la 
progenitors, thus it takes ~ 10—100 Myr for the CO WD 
to accrete up to the Chandrasekhar mass. Therefore, the 
delay time is: tevoi + (~ 10 — 100) Myr (-1- Iq^ in the 
case of helium WD donors)= 0.2 — 2 Gyr for a typical 
AM CVn SN la progenitor to produce a SN la. 

The SDS rate maintains a nearly constant SN la rate 
through about 6 — 7 Gyr and then drops by almost an 
order of magnitude. The increased rate at shorter de- 
lay times is due mainly to systems with evolved donors 
which encounter, on average, RLOF on shorter (evolu- 
tionary) timescales than their main sequence counter- 
parts for which the RLOF is encountered on longer (mag- 

^ a ~ 0.3 — IR0, where the primary is a CO WD and the com- 
panion is either a hehum star (A/hc ~ 0.4 M©) or a helium WD 
(AfeoWD ~ 0.3 Mq). 
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netic braking) timescales. 

The resulting rate of potential DDS SNe la varies sub- 
stantially with time. At early times (i ^ 1 Gyr) the 
rates are very high 0.01 yr^^, and then they gradually 
decrease to reach ^ 0.0003 yr^^ at late times {t > 10 
Gyr). The observed rates for elliptical galaxies are esti- 
mated at the level of 7^obs 0. 0018 ± 0.0006 vr~^ pe r 
unit (10^°L|) of blue luminosity ( [Cappellaro et al.|l999 |. 
As the blue luminosity of elliptical galaxies declmes witii 
time (after an early star formation episode), the rates 
presented in the top panel of Figure 2 should be cor- 
rected downwards at early times, while at later times 
they should be increased if our rates are to be compared 
with those of typical ellipticals. Obviously, the burst of 
star formation on the order of 6 x 10^*^ Mq would produce 
a blue luminosity larger than 10^"Xq, while 10-15 Gyr 
after the episode when stars more massive than ~ 1 Mq 
have formed remnants and are no longer contributing to 
the galaxy's light, the blue luminosity is smaller than the 
normalising value. Since we do not really know the dis- 
tribution of age of the galaxies in the observed sample of 
ellipticals that were used in the SN la rate estimate, we 
do not attempt to correct our synthetic rates for the evo- 
lution of blue luminosity and we do not compare them 
directly to the observed rates of Cappellaro et al. ( 1999 1. 
However, we note that the observed rate is consistent 
with our predicted rates for the DDS progenitor, while 
the predicted rates for other progenitors (SDS and AM 
CVn) seem to be significantly too low. 

In Figure 2 (bottom panel) , we show the SN la rates for 
the spiral galaxy model of Model 1. It is found that DDS 
rates of SNe la at the current epoch are 0.002 yr"-"^. At 
first, the DDS rate increases with time (after the onset 
of star formation), then remains approximately constant 
until the star formation stops leading to an overall de- 
cline in the rate. This behavior reflects the specific shape 
of the delay time distribution for the DDS combined with 
the SFR for our spiral galaxy model. The rates for SDS 
and AM CVn progenitors are much smaller and at the 
level of ~ 10~^ yr~^. SDS progenitors can generate SNe 
la long after star formation has ceased (long delay times) , 
while AM CVn events disappear shortly after the star 
formation has stopped (short delay times). For compar- 
ison, over-plotted are e mpirical rates of SNe la. The 
rates were adopted from [Cappellaro et al.| ([1999') for a 
Milky Way type spiral (sbc-Sd) with a blue luminosity 
of 2 X 10^° Lq, and the rates are 7^obs = 0.004 ± 0.002 
SN la yr~^. The DDS rate alone is consistent with the 
empirical rate of SNe la. The SDS and AM CVn SN la 
rates do not even come close to the empirical rate, and 
their addition to the DDS rate does not significantly af- 
fect the overall rates at any epoch. We note that our 
mass normalization which implies a constant star forma- 
tion history for 10 Gyr results in a SFR at the level of 
6 M0 yr-i. The global SFR in the MW may be some- 



up, being consistent with the DDS scenario as the ma- 
jor SN la contributor in MW-li ke spiral galaxies, as long 
as the Cappellaro et al. (1999) rates are being used for 



what lower: ^ 3.5 M 



yr 



(Cox 



2000 



O'Shaughnessy 



et al. 2008), and in that case the DDiS rates are only 



margi 



nally consistent with the observed, Cappellaro et al. 



([1999| rates. On the other hand, it has been suggested 
that the SFR of the MW has been decreasing with time, 
only reaching ~ 3.5 M^ yr~^ at the current epoch (Nele- 
mans et al.|2001 2004 see sect. 2.2). If such an estimate 
had been used the average SFR of the MW is found at 
the level of ^ 8 Mq yr^^, and our results would scale 



comparison. 

Model 2. There is a marked decrease, by nearly a fac- 
tor of 2, in the total number of SNe la progenitors in our 
model with decreased CE removal efficiency. The overall 
decrease is due to the fact that the most dominant poten- 
tial channel, the DDS, is only ^ 50 % as efficient, since a 
larger fraction of binaries will merge in the common enve- 
lope phase rather than surviving the CE to subsequently 
form a double white dwarf. 

In Figure 3 (top panel), SN la rates are shown for our 
Model 2 elliptical galaxy (acE x A = 0.5, instantaneous 
starburst at t = 0). DDS SNe la progenitors continue to 
outnumber the SDS and AM CVn progenitors, however 
there are some notable differences. For short delay times 
t ^ 1 Gyr, the DDS rates are nearly a factor of 2 lower 
than they are for Model 1. Then at later times, the Model 

2 DDS rates are a factor of ~ 3 below those of Model 1, 
reaching ^ 9 x 10~^ yr""'^ at delay times of 10 Gyr (vs. 

3 X 10"^ yr-i for Model 1). Despite the lower DDS rates 
of Model 2, potential progenitors are found at all delay 
times, as they are in Model 1. The CE efficiency of Model 

1 allows for DDS progenitors to be drawn from a wider 
range (going to smaller values) among the distribution of 
initial separations, where as progenitors with small ini- 
tial separation in Model 2 are removed from the DDS 
population in mergers during the CE phase. The Model 

2 SDS channel is more efficient (by a factor of 3) than the 
Model 1 SDS channel since post-CE WD + MS binaries 
are found on closer orbits. One major difference between 
the elliptical galaxy la rates of Model 1 and Model 2 is 
that the SDS rates match those of the DDS rates for delay 
times ~ 2.5 — 5.5 Gyr (~ 0.0002 yr~^). Since the major- 
ity of the donors are evolved stars (giants or sub-giants) , 
the delay involves two components: the main sequence 
lifetime of the donor (a few Gyr for a donor to become 
a giant) and the accretion timescale, over which the pri- 
mary WD can increase its mass to the Chandrasekhar 
mass (10-100 Myr). Thus the main sequence lifetime of 
the donor is what sets the delay times for the SDS DTD. 
Since the stars are found on closer orbits after the com- 
mon envelope, RLOF is encountered between the WD 
and the non-degenerate companion more often in Model 
2 (typically when the donor is a sub-giant). Rates of 
potential AM CVn SN la are lower than those of Model 
1 for elliptical galaxies and are at the level of ~ 0.0002 
yr~^. For the majority of the progenitors the delay times 
are very short, so as in Model 1, these type of events are 
expected only in young host galaxies or in regions with 
ongoing star formation. 'Fast' AM CVn progenitors are 
more rare in Model 2 since these systems more readily 
merge in one of the two common envelope phases that 
lead to the formation of these progenitors. In contrast to 
Model 1, there is a small contribution of the 'slow' AM 
CVn progenitors (long delay times) in Model 2. The AM 
CVn channel is outnumbered by both the SDS and DDS 
channels at all epochs in the Model 2 elliptical galaxy. 

In Figure 3 (bottom panel), SN la rates are shown 
for our Model 2 spiral galaxy. It is found that DDS 
rate of SNe la at the current epoch (10 Gyr) is 0.001 
yr~^; a factor of two below that of the Model 1 spiral 
galaxy. At first, the DDS rate increases and then re- 
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mains fairly constant until star formation ceases at 10 
Gyr. At all epochs (particularly during star formation), 
the DDS channel rates significantly outnumber the SDS 
and AM CVn channels, but to a lesser degree than when 
compared with Model 1. The SDS channel also exhibits 
a relatively constant rate at times later than ~ 1 Gyr. At 
the current epoch the SDS rates are at the level of 0.0002 
yr~^; a factor of ~ 5 below the DDS rates. The rates 
arising from AM CVn progenitors are fairly negligible 
at all epochs (< 10"'* yr"^). Even when the rates from 
all three progenitor channels in the considered galaxy 
model are combined (0.0012 SNe la yr~^), the SN la rate 
at the cu rrent epoch falls below t he empirically-derived 
rate from Cappellaro et al. ( 1999 1 by roughly a factor of 
two. 

It is worth com paring our mod e l gala xy rates to the 
rates presented in Mannucci et al. ( 2005 1 . In that study, 
Mannucci et al. (2005p derive SN rates for galaxies of 
various morphological types, and present the SN rates 
in SNuM (as inferred from K-hand luminosity measure- 
ments) . We cannot directly compare our ra tes in Tables 
1 & 2 to those of Mannucci et al. ( 2005 1 since we do 
not know the exact ages of the galaxies in their sam- 
ple. However, we note that our Model 1 elliptical galaxy 
SN la rate at 10 Gyr^ is 0.005 SNuM (see Table 1 DDS 
rates) , which is nearly a factor o f 10 below the SN la rat e 
in E/SO galaxies presented by Mannucci et al. (20051: 



0.044;[;;°i| SNuM (see their Table 2). We note as well 
that for the same model galaxy, at t = 500 Myr (shortly 
after a burst of star formation at t = 0) we obtain a SN 
la rate of ~ 0.18 SNuM (mostly via the DDS channel 
with some contribution from AM CVn), which is about 
a factor of 2 lower than the range of SN la rates found for 
star forming Irre gular galaxies in [Mannucci et al. ( 2005 1 
The 



(0.771°:^?). 



Mannucci et al. (|2005|) SN la rate for 
SNuM, 



SOa/b spirals is found to be 0.065^^^ SNuM, which 
matches our Model 1 spiral galaxy rates at 5 Gyr (0.065 
SNuM; mostly DDS & AM C Vn progenitors) The SN 
la rate of Sbc/d spirals from Mannucci et al. ( 2005[ ) is 
0.17tg;o63, and matches our Model 1 spiral SN la rate 
only at very early times (~ 0.5 — 1 Gyr, see Table 1). 

For the Model 2 elliptical galaxy, we find a low SN 
la rate of 0.001 SNuM at 10 Gyr, which is over an 
order of magnitude be low the rate for E/SO galaxies in 



014 SNuM). The rate at 



Mannucci et al.| ( |2005[ ) (0.044^° °^^ 
500 Myr for the same galaxy is 0.14 SNuM (see Ta- 
ble 2); mostly arising from DDS progenitors with some 
contribution from AM CVn and SDS. This rate is a fac- 
tor of a few belo \ y the la rates for Irregular galaxies in 



Mannucci et al. (2005) (0.77I^;3i). In comparing our 
Model 2 spiral rates, we find that our SN la rate at 3 Gyr 
(~ 0.07 SNuM; mostly DDS progenitors with some con- 
tr ibution from SDS) is wi thin the range of rates presented 



in Mannu cci et al.| ( |2005 | for SOa/b spirals (0.065;^ |525)' 
while only our spiral rates for Model 2 at very early times 
(< 1 Gy r; 0.13 SNuM, DPS) a re high enough to match 
those of 'Mannucci et al. (20051 for Sbc/d type galaxies 



(0.17j;o!og|). We note that in general, our rates (per unit 



mass) are lower than those of Mannucci et al. (20051, in 



^ A typical age for local ellipticals; see Mannucci et al. 2005 
Section 7. 



dicating that perhaps other channels leading to the for- 
mation of SNe la should be considered in evolutionary 
studies (e.g., single stars, sub-Chandrasekhar mass SNe 
la). We also note however that our predicted rates as 
a function of time (at least for the D DS) are consistent 
with t he observed rates presented by [Cappellaro et aL\ 
( [1999| . 

We find that in general the DDS outnumber SDS and 
AM CVn progenitors. This effect is somewhat more pro- 
nounced in Model 1 (a factor of > 10) than in Model 2 
(a factor of ~ 5), but the reason why is clear for both 
models. The occurrence rate of a CO-CO WD merger 
with a total mass > 1.4 Mq (DDS) is higher than that 
of building up a CO WD's mass to ^ 1.4 Mq via stable 
mass transfer (SDS) in a binary. Formation efficiencies 
in both of the above cases are very low, after all SNe la 
are rather rare events. However, relatively speaking it is 
easier to find a pair of two CO WDs (DDS), each with a 
mass of > 0.7 Mq which is a typical mass for CO WDs, 
than it is for a CO WD to double its mass through accre- 
tion (SDS/ AM CVn). This finding is a consequence of 
the recent updates on accumulation physics calculations 



al. 2b07|) that we have adopted in our evolutionary study. 
caIR 



(jHachisu et al. 1999 Kato & Hachisu 



1999 Nomoto et 



BasicaUy, the accumulation onto a WD is hampered by 
a number of processes that tend to remove matter which 
is transferred from the companion in a close binary (e.g., 
nova explosions. He-shell flashes, optically thick winds 
from the surface of an accreting WD), in some cases 
leading to the disruption of an accreting WD before it 
reaches the limiting Chandrasekhar mass; i.e., accretion 
of ^ 0.1 M0 of a He- rich layer which ignites and disrupts 
the underlying WD (sub-Ch andrasekhar mass SN la, see 



e.g., Kato & Hachisu 19991. Since, (i) there is a rather 



narrow range of mass transfer rates which may lead to 
efficient accumulation and (ii) there are not that many 
binary configurations (and we have considered the entire 
range for our adopted evolutionary model) that can sus- 
tain mass transfer for a prolonged period of time, the 
SDS and AM CVn channels are found to produce SNe la 
at very low rates. 

5. COMPARISON WITH OTHER STUDIES 



The recent theoretical study of Hachisu et al. ( 2008 1 
finds a SDS delay tim e distribution wh ich follows a power 
law 



In their study, 



Hachisu et al] (20(18 

a new mass stripping ett'ect (based on |H 

(1999l), where in the case of high mass transfer rates 



incorporate 
" PSachisu et al 



the WD blows an optically thick wind strong enough to 
'strip' material from a main sequence or giant donor. 
This effect in return stabilizes mass transfer, enabling 
the binary to avoid a CE phase even in the case of a rel- 
atively massive 6 Mq) donor. The result is that the 
WD can accrete stably up to the Chandrasekhar mass, 
with a wider range of poten tial progenitor donor ZAMS 
masses: 0.9 - 6 Mq in Hachisu et al.j (|2008| vs. 0.7 - 2.7 
Mq in our current study. Even though we allow for SN 
la progenitors to form from any initial mass spanning 
the initial mass function, our SDS SNe la only derive 
from binaries involving low-mass donors since we do not 
take into ac count this stripping effect. We note that the 
Hachisu et al. (2008) model predicts the presence of a 
thick disc of hydrogen-rich circumstellar material around 
the SN la progenitor. If such a circumstellar torus were 
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present around the majority of progenitors of SNe la, one 
would likely expect to observe hydrogen in their spectra, 
though to date less than 1% of SNe la h ave shown any 



signature of as sociated hydrogen (see e.g., Han & Podsi- 
[adlowski 2006 



|Han fc Podsiadlowski (2004) investigated SNe la pro- 
genitors from WD + MS and WD + evolved binaries 
with very specific binary configurations using population 
synthesis. In that work, they do not present rates of SNe 
la derived from other possible formation channels of SNe 
la so we cannot compare DDS rates. The Galactic rate of 
SDS SNe la for their model which most closely matches 
our standard (Model 1) parameters is 6 x 10~* yr~^. 
This value is an order of magnitude above our SDS spiral 
galaxy Model 1 rate of ^ 6 x 10~^ yr~^, though is still 
nearly an order of ma gnitude lower than the empirical 



Cappellaro et al. ( 1999 4 x 10"^ 



SN la Galactic rate of 

yr~^). The rate for the |Han &: Podsiadlowski 



model which most closely matches Model 2 is ~ 10 



(20041 

la yr~^ , which is close (though still below) the Galactic 
rates of Cappellaro et al. ( 1999 1, and is a factor of 5 times 
higher than our Model 2 SDSTrates (^ 2 x 10""^ yr ~^). 



It was pointed out by Han & Podsiadlowski ( 2004 ) that 
their prescription for hydrogen accumulatio n is more effi- 
cient than that used by other authors (i.e., 'Yungelson fc] 
Livio|1998[ ). It is also more efficient than the prescription 



we have adopted in this study. The range of hydrogen 
accretion rates onto WDs which leads to stable burning 
(see section 2) and efficient mass accumulation is uncer- 
tain, and it is possible that stable hydrogen burning may 
occur for a wider range of accretion rates, in turn allow- 
ing for higher S Ne la rates following from the SDS chan- 
nel as allowed in Han & Podsiadlowski (20041. However, 
comparison of model hydrogen-accreting WDs on the H- 
R dia gram with supersoft X-ray sources (Nomoto et al 



|200 

bur 



2007| indicates that the prescription for stable hydrogen 



urning in a thin shell (and adopted here) is consistent 
with observations. 

It is worth no ting that the delay time distributions 
of|Greggio (|2005[ ), derived using analytical formulations, 
produce a DDS DTD shape which is similar to ours: 
peaked at short (< 1 Gyr) delay times, followed by a 
smooth drop-off as a function of time, due to the de- 
pendence of the delay time on the timescale associated 
with gravitational wave emission. A similar trend is also 



found for the delay times of DDS progenitors in Yungel- 



son fc Livio| pOOO| Fig. 2). The jGreggiol ( |2005p stud~ 
also determined that the shape of the DI D arismg from 
SDS progenitors was more flat than when compared to 
that of the DDS, and that the SDS delay time depended 
upon the main sequence lifetime of the secondary star 
(see their section 5), which is consistent with our find- 
ings. 

Delay tim es of SNe la were calculated by Belczynsiol 
et al.| ( 2005 1 . It was found that the merger of two white 



dwarfs was consist ent with an einpiri cal delay time esti- 
mate of - 3 Gyr (Strol ger et"aD '2004)^, and that WDs 
accreting from non-degenerate stars could potentially ex- 
plain the observed delay times if a low common envelope 
efficiency is used (acE x A = 0.3) and if it is presumed 
a priori that WD mergers contribute negligibly or not 

Delay times were compu ted by adopting a co smic star forma- 
tion history based on that of |Madau et al.| | 1998| . 



at all to the SNe la population. The above results were 
obtained with an earlier version of the StarTrack code; 
the code was recently up dated to include the m ost recent 
accumulation rates (e.g., Nomoto et al.|[2007 l. All other 
recent revisions relevant for low- and intermediate-mass 
binary evolution in the code are described i n detail in 
Belczynski et al. ( 2008 1 . Additionally, in the Belczynski 
( |2005p study, very different criteria were adopted 



et al. 



for sNe la; the majority of their supernovae originated 
from sub-Chandrasekhar mass events and it was permit- 
ted that the merger of two WDs of any type (including 
ONeMg WDs) with a total mass exceeding 1.4 Mq led 
to a SN la. In the current study we also note that the 
number of sub-Chandrasekhar mass events exceed the 
number of Chandrasekhar mass progenitors. The weak- 
ness of the sub-Chandrasekhar model is that much of the 
outer (and fastest moving) material is believed to b urn to 
nickel with ver y few intermediate mass elements (Livne 
fc Arnett|1995 ). The resulting spectra do not match cur- 
rent observations of normal supernovae, and the best fits 
are for sub-luminous supernovae. Only ~ 6 sub-luminous 
SNe la were recen tly reporte d (Kasliwal et al. 20 0^ as 
comp ared with 36 ( Riess et al.|1998| ) or 42 (Strolger et al. 
2004 1 normal ones discovered only in the Hubble surveys. 



6. DISCUSSION 

We have evolved single and binary stars using the pop- 
ulation synthesis code StarTrack, and have analyzed 
the resulting delay times (time from binary formation 
at t = to SN la) of potential Type la supernovae. 
We have considered possible SNe la progenitors aris- 
ing from three formation channels: Double Degenerate 
Scenario (white dwarf mergers). Single Degenerate Sce- 
nario (hydrogen-rich accretion on to a WD) and the AM 
Canum Venaticorum scenario (helium-rich accretion on 
to a WD). Additionally, we have computed SN la rates 
for two galaxy types: an elliptical galaxy with a starburst 
at t — and a spiral galaxy with a constant SFR, and 
in each case we have tested the impact of the common 
envelope removal efficiency for two different parameter- 
izations: acE x A = 1 (Model 1), and acE x A = 0.5 
(Model 2). 

Our SN la rates (century)^-'^ (lO^'^Mo)^^ have been 
presented in Tables 1 and 2. We reiterate that the rates 
which we have derived in this work are 'local' (no red- 
shift evolution) rates. It is still interesting to note that 
in calculating the volumet ric SN la rate out to z = 0.12 
fitting the A -I- B model, [Dilday et al] ( |2008[ ) found A 
= (2.8 ± 1.2) X IQ-^'^ SNe la yr'^ M~\ which is most 
similarly matched by our DDS Model 1 rates for Milky 
Way-like spirals (- 3 x lO^^^ SNe la yr"! M~^ at 10 



Gyr). This 2-component or "A-|-B" model (see Mannucci 



et al.| 2006 Scannapieco fc Bildsten 2005) assumes that 
the SN la rate is a function of stellar mass density and 
the SFR, and allows for fitting prompt and tardy SN la 
populations. We note that for Model 1, DDS progenitors 
dominate the overall SN la rates, whether it is for ellipti- 
cal or spiral hosts. For Model 2, the overall rates of SNe 
la are a factor of two lower than in Model 1, though SNe 
la progenitors arising from the DDS and SDS channels 
are found in equal numbers for delay times 2.5 — 5.5 Gyr. 
For all models considered, we expect both very short 
1 Gyr) and long (- 10 - 15 Gyr) SN la delay times if 
they originate from DDS progenitors. Could such a pop- 
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ulation explain the bimodal distribution of delay times 
derived from some observations? In principle, one may 
expect such a result; if DDS progenitors are dominant 
in both old and young galaxies, the empirically-derived 
delay times may appear to be bimodal and erroneously 
point toward two different progenitor populations. How- 



ever, the reported bimodal DTD among SNe la (Man 



nucci et al. 



20061 i-e., 50 % of SNe la having 'proSpT 



delay times < 100 Myr) is not reproduced in any of our 
models. For our standard model, only ~ 5% of systems 
have delay times below 100 Myr, with 50 % of our SNe 
la occurring within ^ 800 Myr since the starburst. The 
model which comes closest to reproducing such a bimodal 
delay time distribution is the DDS channel of Model 2, 
where there is a large fraction of SNe la progenitors with 
short delay times (50% with t < 500 Myr). If in fact 
such a large fraction of SNe la occur within 100 Myr 
of star formation, it may pose a very interesting prob- 
lem for binary evolution to explain. Some alternatives 
are already being cons idered, for example a single star 
SN la progenitor (Iben fc Tutukov] |1984] [Maoz 2008| . 
To properly approach this issue one needs to fold our 
evolutionary calculations with the cosmic star formation 
history and distribution of galaxy types and mass as a 
function of redshift. However, it was pointed out that 
the constraints on observational DTDs which incorporate 
convolution with assumed cosmic star formation history 



may not be as strong as they are claimed to be (|F6rster 
eral][2006| 



Tot our standard (Model 1) set of evolutionary param- 
eters, the DDS dominates (roughly 90%) the rate of Type 
la SNe in both spiral and elliptical galaxies. The rate of 
DDS SNe la is consistent with the observed SN la rate 



of |Cappellaro et al.| ( [1999| . The DTD of potential DDS 
SNe la follows a smooth power law distribution through 
a Hubble time. The SDS DTD is mostly flat throughout 
a Hubble time, while the AM CVn SNe contribute events 
mostly at short delay times, but neither contribute signif- 
icantly to the total SN la rate for spiral-like galaxies. The 
DDS DTD, following a power la w, is consistent with the 



findings of Totani et aL] ( 2008 1 , whose observationally- 
derived delay time of SNe la in intermediate-redshift el- 
liptical galaxies follows a featureless power law oc for 
0.1 < t < 10 Gyr. 

For Model 2, the DDS dominates (83%) the rate of 
Type la SNe in both spiral and elliptical ga laxies. The 
DBS rate is below the observed SN la rate of |Cappellaro| 
et al.Ml999] ) by a factor of two. The DTD of potential 
DDS sNe la does not exhibit as strong of a power-law 
shape as Model 1, since relatively more white dwarf bi- 
naries merge at early delay times, and more would-be 
DDS progenitors merge during the CE phase and thus 
never produce a SN la. The SDS DTD is most promi- 
nent ~ 2 — 6 Gyr (at which times the SDS rates match 
those of the DDS) and at later times the SDS contributes 
very little to the overall rates. The AM CVn SNe con- 
tribute some events mostly at short delay times, with the 
rates being at a very low level at long delay times. 

As noted previously, all DDS SNe la progenitors are 
the result of a merger of a CO-CO WD binary. Mergers 
between CO-He WDs are relatively less common than 
He-He or CO-CO WD mergers since often times mass 
transfer is stable, and upon reaching contact the CO-He 
binary will enter an AM CVn phase rather than coalesce. 



In any case, mergers of He-He or CO-He WDs are not 
massive enough to lead to SNe la, though they may lead 
to other very interesting phenomena, as they are likely 
precur sors to he lium-burning hot subdwarfs or R CrB 
stars (Webbink 1984). There is a marked decrease in 
the number of double white dwarfs which are formed in 
Model 2 relative to Model 1 due to a heightened number 
of mergers which are encountered during the CE phase. 

Typically, population synthesis calculations produce 
the right number of Chandrasekhar mass WD- WD merg- 
ers, yet usually produce an order of magnitude too few 
SDS supernovae to explain the observed rate estimates 
(see LiviolpOOT for a review). As it turns out, most hy- 
drodynamical white dwarf merger calculations result in 
the less massive white dwarf being disrupted in a few 
orbits, causing the more m assive white dwarf to accrete 



at an extrem ely high rate (iRasio & Shapiro 1995 Benz 
ligh 



et al.| 19891) . At very high accretion rates, the white 
dwarf mergers are believed to collapse to form a neu- 
tron star, not produce a SN la thermo nuclear explosion 



jMochkovitch fc Livio|| 1990l |Saio fcTTomoto 1985 , 1998' 
Woosley fc Wea ver 199^~S ome new studies have found 
that some mergers might produce Type la super novae, 
but the mer ger conditions must be severely refined (Yoon 



et al][2 007l. 
!Eh 



Either population synthesis studies with the adopted 
accretion physics are missing active progenitor paths for 
Type la supernovae, or the physics of merger calcula- 
tions are incorrect. All of our results and conclusions 
were based on one evolutionary model with a specific 
(albeit the most updated) set of accumulation rates and 
have been presented from the standpoint of population 
synthesis (i.e., thus far we have i gnored the inform ation 
provided by merger calculations). Yoon et al. (20071 pre- 
dicted that only DDS systems with small mass ratios 
g < 0.4 can produce SN la explosions, while the rest will 
end up in accretion induced collapse and neutron star 
formation. If we had adopted this as an additional cri- 
terion the implications would have been rather dramatic 
for our results. At the time of the merger, only ~ 0.2% 
(0 %) of the DDS systems from Model 1 (Model 2) have 
q < 0.4, since a low q usually leads to stable mass transfer 
and not a merger. If the DDS rates had been decreased 
so drastically, none of our calculated SN la progenitors 
would be able to match the observed rates. 

If we had relaxed our assumption on Chan- 
drasekhar mass explosions and we had incorporated sub- 
Chandrasekhar mass explosion models - ignition of a de- 
generate layer of H e-rich material accumulated on the 



surface of a WD^ ([Woosley fc Weaver 
Hachisu 1999 1 - the Galactic rates for AiV 
crease from 



1994 



Kato fc 



CVn would in- 

_ 10"^ yr-i to ~ 10-3 yr-i for Model 1 and 
to ~ 5 X 10"^ for Model 2. We note that in this estimate 
we have allowed sub-Chandrasekhar mass explosions to 
occur only for white dwarfs with masses M^d > 1 M0, 
since for lower WD masses the explosions would not look 
like those observed for ty pical Type la supernovae (e.g., 
Hillebrandt fc Niemeyer||2000) . 

If indeed white dwarf mergers (DDS) cannot lead to 
thermonuclear explosions, one needs to consider alter- 
natives for increasing the rates of SDS and AM CVn 

* The sub-Chandrasekhar mass model falls under our AM CVn 
formation channel. 
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progenitors, cither by moving away from the standard 
evolutionary model, or widening the range for efficient 
accumulation rates onto white dwarfs. 
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TABLE 1 

Rates of SNe Ia Progenitors (SNuM) for Model 1 



Rate [(100 yr)-i 10^°M~^] as a Function of Time. 



DDS 


Elliptical 


Spiral 


0.5 Gyr 


1.6 X 10-1 


2.0 X 10-1 


3 Gyr 


2.3 X 10-2 


8.0 X 10-2 


5 Gyr 


1.2 X 10-2 


6.0 X 10-2 


10 Gyr 


~ 5 X 10-^ 


3.3 X 10-2 


SDS 






0.5 Gyr 


< 10-3 





3 Gyr 


^'s X 10-3 


< 10"^ 


5 Gyr 


~ 1 X 10-3 


~ 10-3 


10 Gyr 


< 10-3 


~ 10-3 


AM CVn 






0.5 Gyr 


2.2 X 10-2 


~ 10-2 


3 Gyr 


< 10-3 


~ 5 X 10-3 


5 Gyr 


< 10-4 


~ 4 X 10-3 


10 Gyr 


~ 


~ 1 X 10-3 



TABLE 2 

Rates of SNe Ia Progenitors (SNuM) for Model 2 



Rate [(100 yr)-i IQIOM"!] as a Function of Time. 



SDS 





Elliptical 


Spiral 


0.5 Gyr 


1.3 X 10-1 


1.3 X 10-1 


3 Gyr 


~ 5 X 10-3 


6.0 X 10-2 


5 Gyr 


~ 3 X 10-3 


3.5 X 10-2 


10 Gyr 


~ 10-3 


1.4 X 10-2 


0.5 Gyr 


~ 2 X 10-3 





3 Gyr 


~ 5 X 10-3 


~ 8 X 10-3 


5 Gyr 


~ 3 X 10-3 


~ 6 X 10-3 


10 Gyr 


~ 10-4 


~ 3 X 10-3 


CVn 






0.5 Gyr 


~ 5 X 10-3 


< 10-3 


3 Gyr 





~ 10-3 


5 Gyr 





~ 10-3 


10 Gyr 


< 10-4 


< 10-3 




Fig. 1. — Delay time distribution showing relative contributions from the three SN la formation channels considered in this work for 
our elliptical galaxy (instantaneous burst of star formation at t = 0): DDS (blue dashed), AM CVn (green dot-dash) and SDS (red solid). 
Average and median delay times are indicated. Top panel: Model 1, acE x A = 1.0. The thin dotted line represents a power law function 
of form f{t) = 100 and follows the DDS delay time distribution reasonably well. Bottom panel: Model 2, acE x A = 0.5. The power 
law is shown for comparison. 
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Fig. 2. — Rates of Type la supernovae (number per year) for the Model 1 galaxies. The DDS (blue dashed), AM CVn (green dot-dash) 
and SDS (red solid) channel rates are shown. Top panel: elliptical galaxy whose total mass in formed stars is Mtot = 6 X 10^*^ Mq 
(delta function starburst at t = 0). Bottom panel: spiral (constant star formation) galaxy whose total mass in formed stars at 10 Gyr is 
Mtot = 6 X 10^'^ Mq. Approximating the star formation history of the MW to be constant, the current Galactic SN la rate (shown with 
vertical dotted line) considering the DDS alone is 0.002 yr~^. This is consistent with the empirical rate as indicated by the shaded region 
in the plot (0.004 ±0.002 yr- MCappellaro et al.|l999| . 
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Fig. 3. — Same as Figure 2 for the Model 2 stellar population. Top panel; elliptical galaxy. Bottom panel; spiral galaxy. Approximating 
the star formation history of the MW to be constant, the current Galactic SN la rate considering the DDS alone is 0.0009 yr~^. The 
combined (DDS + SDS + AM CVn) rate is 0.001 yr'^, which is below the empirical rate estimate of Cappellaro et al. (1999; 0.004 ±0.002 



